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Abstract 

We numerically investigate multiphoton ionization of excited hydrogen atoms 

in magnetic and microwave fields when up to Nj = 600 photons are required 

for ionization. The analytical estimates for the quantum localization length 

in the classically chaotic regime are in agreement with numerical data. The 

excitation is much stronger as compared to the case with microwave field only 

due to the chaotic structure of eigenstates in magnetic field. 
P.A.C.S.: 36.10.-k, 05.45.+b 
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In 1982 Shushkov and Flambaum discussed the effect of weak interaction enhancement 
due to a complex structure of ergodic eigenfunctions in nuclei. The basic idea of this effect is 
that in complex systems an eigenfunction, represented in some basis, has a large number M 
of randomly fluctuating components so that their typical value is 1/ \/M. Due to this, the 
matrix elements for interparticle interaction are Vint ~ 1/ a/M while the distance between 
mixed levels is AE ~ 1/M. As a result, according to the perturbation theory, the admixture 
factor ?7 is strongly enhanced: "q ~ Vint/^E ~ v^M as compared to the case in which 
eigenfunctions have only few components (M ~ 1). This effect was investigated and well 
confirmed in experiments with weak interaction enhancement for scattering of polarized 
neutrons on nuclei [0. Recently a similar effect of interparticle interaction enhancement 
was discussed for two interacting particles in disordered solid state systems 0. Here, a 
short range interaction produces a strong enhancement of the localization length leading to 
a qualitative change of physical properties. This shows that the effect is quite general and 
can take place in different systems. 

In this Letter we investigate the possibility of similar enhancement in atomic physics for 
atoms interacting with electromagnetic fields. Such process becomes especially interesting 
for highly excited atoms (hydrogen or Rydberg atoms) in microwave fields where absorption 
of many photons is necessary in order to ionize electrons. Until now this problem was studied 
only in the case in which the electron dynamics, in absence of microwave field, is integrable 
0. In this case strong ionization is possible due to onset of chaos at sufficiently strong 
field intensity. As it is known, above the classical chaos border ionization proceeds in a 
diffusive way and quantum interference effects can lead to localization of this diffusion and 
thus suppress ionization. 

A quite different situation, which was never studied neither numerically nor experimen- 
tally, appears when the electron's motion in the atom is already chaotic in the absence of 
microwave field. An interesting example of such situation is an hydrogen atom in a strong 
static magnetic field. The properties of such atoms have been extensively studied in the 
last decade and it has been shown that the eigenfunctions are chaotic, and that several 
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properties of the system can be described by Random Matrix Theory. Due to that one can 
expect that the interaction of such an atom with a microwave field will be strongly enhanced 
so that the localization length will become much larger than the corresponding one in the 
absence of magnetic fields. As a result, the quantum delocalization border, which determines 
the ionization threshold, will be strongly decreased. 

The investigation of classical dynamics and some preliminary estimates for the quantum 
localization length in such a case have been given in a recent paper Here we discuss 
the quantum dynamics and present the results of numerical simulations which confirm the 
chaotic enhancement of quantum excitation as compared to the case without magnetic field. 
Our studies also show a number of interesting features which arise in this model in the 
adiabatic regime when the microwave frequency is much smaller than the Kepler frequency. 

We consider the case in which the electric and magnetic fields are parallel. In this case 
the magnetic quantum number m is an exact integral of motion and here we set m = 0. The 
Hamiltonian writes 



where ul = B/c = B(T)/Bq is the cyclotron frequency, Bq = 2.35 x lO^T, e and u are the 
field strength and frequency respectively (atomic units are used). As it is known @J^, in 
the absence of microwave field, the classical motion becomes chaotic for uj^riQ > 1 and no 
visible islands of stability are present for uinl ~ 9. For cul^o ~ 3 some islands of stability 
exist but their size is small. 

The turn on of microwave field leads to diffusive energy growth with classical diffusion 
rate per unit time Db = (AE)^ / At. The dependence of Db on parameters e, u has been 
found in §]: Db/Dq ^ Xi^o i^o < 1), Db/Dq ^ X2/t^o^^ (^^o > 1), where ujq = unl, 
Dq = is the diffusion rate in the chaotic regime for i? = and ujq = 1 and Xi? X2 

are two constants, weakly dependent on magnetic field (numerically, xi ~ 18; X2 ~ 2 at 
^L^l = 9.2 and Xi ~ 25, X2 ~ 1 at uj^nl = 3). The above estimates for Db give the 
asymptotic behavior of the diffusion rate for very small and very large uq, while the actual 



values of Db were determined from numerical simulations of the classical problem. 

In the quantum case the interference effects can lead to localization of this diffusion 
PJ^. The localization length in number of photons is is proportional to the one-photon 
transition rate T and to the density of states ps coupled by these transitions 0: is ~ Tps- 
The transition rate F can be derived from the classical diffusion rate: F ^ Db/oj"^- We recall 
that for the case _B = the localization length (i^ aX ujq = 1 is = S.SegnQ ~ Dqti^pq [0, 
where po = ""-o is the density of states and eo = en^. Due to Coulomb degeneracy and to 
the existence of an additional approximate integral of motion the density po is Uq times 
smaller than the number of levels in one unit energy interval. As a result 

„ „ Db Pb 
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where it is assumed £b > 1 and upB > 1. According to our quantum data ps/Po ~ ^o/('^l'^o) 
for UluI > 1. More exactly ps/po = 0.34^0 (for ^L^l = 3, Uq = 60) and Pb/Po = 0.14no (for 
f^L^o — 9-2; ^0 = 60). The dependence pb ~ I/^^l is due to the oscillatory type behavior 
in p direction in (|l]). The number of photons required for ionization is Nj = nQ/2ujQ and 
therefore for is Nj eigenf unctions are exponentially localized in the number of photons 
N^, namely ipN ~ exp{-\N^\/iB)- 

The value of the localization length £s is strongly enhanced compared to the length 



= 3.3elnl/ujl^^^ at S = and ujq > 1. The enhancement factor (b/^cP^ ~ X2Pb/p> 







X2^o/(i^L^o) ^ 1 is proportional to the initially excited state no- In fact in the presence of 
a magnetic field there is no additional integral of motion [|^ and the number of components 
in the eigenf unctions is increased by a factor M = ps/po- As a result the admixture factor 
rj is also enhanced, namely rf M similarly to the enhancement of localization length in 
disordered solid state models with two particles (^b/^0^ ~ ^ j^y 

The condition Cb = gives the delocalization border above which quantum excitation 
is close to the classical one (both for uq <1 and ujq > 1): 
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For ujq > 1 this value is approximately by the factor (S/rio)^''^ below the delocalization 
border in microwave field only {B = 0) where 

eqo ~ ujI^VVQ^- For uq < 1 the boder is 

eg = {ujL/6.6xiujy^^/no. 

In order to check the above estimates (eqs. (||) and @) we analyzed the quantum 
dynamics following the wave packet evolution in the eigenstate basis at e = 0. Initially 
only one eigenstate is excited with eigenenergy Ex^ Eq = —l/2n^ and uq = 60. In 
our computations we used a total number of eigenstates up to 800 and the evolution was 
followed up to time r = 200 microwave periods. The parameters were varied in the intervals 
0.05 < ujq < 3, 0.002 < Co < 0.02 for oolit-q = 3 and 9.2. For this parameter range, the number 
of photons Nj = no/2u;o required for ionization varies in the interval 10 < Nj < 600. The 
probability distribution f\ over the eigenstates at e = is shown in figs. |l], |^ as a function of 
the number of absorbed photons = {Ex — Eo)/u!. In order to suppress fiuctuations this 
probability was averaged over 10 — 20 microwave periods. For the comparison with classical 
results we also determined the probability Jn in each one-photon interval around integer 
values of A^^. The classical distribution was obtained by solution of Newton equations with 
up to 5 X 10^ classical trajectories and was normalized to one-photon interval. Initially the 
trajectories were distributed microcanonically on the energy surface at energy Eq. 

The typical results in the localized regime are presented in fig. Here the distribution 
reaches its stationary state with a well localized exponential profile. The least square fit 
with /tv ~ exp{—2Nff,/£BN) for > allows to determine the numerical value of the 
localization length £bn which turns out to be in good agreement with theoretical estimate 
(eq. (H)) and is strikingly enhanced compared to the case of zero magnetic field. The plateau 
which appears in fig. |l]b for A^^ > 130 is related to the finite size of the basis and to the 
fact that, according to eq. (|^), the localization length is non homogeneous on high levels 
(is ~ < ~ (Nj - iV^)-"/' for a;o < 1). 

In fig. 1^ the distributions in the delocalized case are shown for ujq = 1 (fig. ^) and 
uo = 0.1 (fig. The delocalization borders in these cases, tq = 0.016 (fig. 0a) and 
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eg = 0.014, (fig. are below the field peak value eo = 0.02. The numerical results show 
a good agreement between classical and quantum distributions in this regime. Even if at 
eo = 0.02 the dynamics starts to be chaotic also at zero magnetic field, the excitation at 
i^L^l = 3 is much stronger due to the chaotic enhancement of electron's interaction with 
the microwave field. For example, the increase of ((AA*",^)^) after 50 microwave periods is 
approximately 55 for the case of fig. |^a while at 5 = it is only 6. 

In order to check the theoretical predictions for the photonic localization length £s we 
analyzed different probability distributions in the localized regime for 0.05 < uq < 3 a.t 
^^L^l = 3 and loltiI = 9.2. The comparison of the numerically obtained lengths isN with 
the theoretical estimate eq. (|[) is presented in fig. ^ Without any fitting parameter the 
theoretical line demonstrates the fairly good agreement with numerical data. Indeed the 
average value of the ratio R = {(-bn /C-b) = 0.81 ± 0.34, where the average is over all values 
iB > 1- The separate averaging over the cases with uq > 1 and uq < 1 gives Ri = 0.98 ±0.30 
and R2 = 0.70 ± 0.26 respectively. 

In spite of the good agreement between theoretical predictions and numerical data we 
would like to stress that a deeper investigation of the problem is required. Especially unusual 
is the regime Uq < 1, which has not been studied up to now from the viewpoint of dynamical 
localization. A number of new questions appear in this regime. For example for cuq > 1 a 
chain of one-photon transitions is clearly seen in the quantum localized distribution (fig. |l]a) 
while for uq < 1 the structure is not visible even though upB > I (fig. |l|b). 

Another interesting question in this adiabatic regime cuq ^ 1 is connected with the 
possibility of analyzing the problem in the instantaneous time basis. In this basis the 
Hamiltonian takes the form H{t) = Ho{t) + dS/dt, where Hq is the Hamiltonian (||) at a 
given moment of time while dS/dt describes the transitions due to the field's variation with 
time (5* is the action of the Hamiltonian (|1|) in which time is considered as a parameter). 
The term dS/dt can be estimated as dS/dt ~ e u sin^ut) 27m^ z eojrt' (see also 0). It 
describes the transitions between instant time levels of the Hamiltonian (|I]), the amplitude 
of which can be estimated as Vg// ~ eun^ / ^/n. The factor ^/n in the denominator appears 
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due to the chaotic structure of eigenstates which leads to a smearing of dS/dt over the n 
states which contribute to the eigenfunctions inside an atomic shell (we assume ujl ~ l/n^). 
Since the distance between levels is 5E ~ l/ps ~ "^o^ it seems that mixing between instant 
levels is possible only if V^// > 5E, giving e^uJoUQ > 1. This adiabatic condition is more 
restrictive than the standard £b > ^ (^o'^o''^ > !)• However the numerical results (fig. ^ 
confirm our estimate (§). For a possible explanation of this discrepancy one may argue 
that the distance between coupled quasi-energy levels is 6E^ u/n and then the condition 
Veff > SE^ gives £b > 1 in agreement with the estimate (@). Another possible reason is 
that in the instant time basis the levels are moving with time and can therefore intersect 
each other giving 6E = 0. 

In conclusion our numerical investigations confirm the theoretical estimates for the pho- 
tonic localization length (0) both for cuq > 1 and ujq < 1. Due to the chaotic structure of 
the eigenstates the quantum delocalization border is strongly lowered compared to the case 
with microwave field only. Since for ujq <^ 1 a. much larger number of photons is required for 
ionization {Nj = 300 in figs. [I|b, |^) experimental observation of localization and verification 
of theoretical predictions should be more easily feasible in laboratory experiments. 
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FIGURES 

FIG. 1. Probability distribution as a function of photon number N^. Quantum distribution 
f\ over the eigenbasis at e = (full line); quantum probability in one-photon interval /tv (circles); 
classical distribution in one~photon interval (dashed line). The straight line shows the fit for the 
exponential decay. 

(a) no = 60, ujq = 1, ujluI = 3, eo = 0.005, Ibn = 3.6, Ib = 3, 180 < r < 200, Db/Dq = 0.49. 

(b) no = 60, Wo = 0.1, oJinl = 3, eo = 0.005, Ibn = 29.5, Ib = 37.5, 180 < r < 200, 
Db/Dq = 0.062. 

FIG. 2. Same as in fig. Q but in the delocalized regime, with eo = 0.02 > e^ 0.015, and 
40 < r < 50. The classical (dashed line) and quantum (circles) distributions in one-photon interval 
are close to each other. 

FIG. 3. The numerically computed localization length (1.bn versus the theoretical estimate Ib 
eq. (§) for w^nQ = 3, cjo > 1 (full circles), uol^^ = 3, cjo < 1 (full triangles), wl'^o = 9-2) '^o ^ 1 
(open circles), uJin^ = 9.2, ujq <1 (open triangles). The error bars obtained from least square fits 
of the localized distributions are also shown. The straight line corresponds to Ibn = ^B- 
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